Suppression, persistence and reentrance of superfluidity in overflowing nuclear systems 
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Based on a microscopic description of superfluidity in overflowing nuclear systems, it is shown that 
continuum coupling plays an important role in the suppression, the persistence and the reentrance 
of pairing. In such systems, the structure of the drip-line nucleus determines the suppression and 
the persistence of superfluidity. The reentrance of pairing with increasing temperature leads to 
additional critical temperatures between the normal and superfiuid phases. 
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Overflowing many-body fermionic systems exist in var- 
ious situations going from the crust of neutron stars 1] to 
ultra-cold atoms [2|, |3J . Interestingly, these systems offer 
the possibility to study the coupling between two fluids 
with very different pairing properties jj, l5| . In such sys- 
tems, one fluid is localized inside an initial container, such 
as for instance a nuclear potential, and a second fluid is 
overflowing towards a larger container. Being in differ- 
ent environments, these two fluids can acquire different 
pairing gaps. In this Letter, we address the question of 
the coupling between the two superfluids and their finite 
temperature properties. 

At the transition between its outer and inner crust, 
neutron stars provide an example of such microscopic 
overfiowing systems, commonly called neutron drip- 
ping [1,] . In the outer crust, nuclei form a Coulomb lattice 
which gets more and more neutron rich as the density in- 
creases. When the maximum number of neutrons that 
a nucleus can sustain is reached, the overproduced neu- 
trons drip out of nuclei. These neutrons populate the 
continuum states and shall be described within the band 
theory (see Refs. [6, 7] and references therein). It should 
be noted that nuclei surrounded by an infinite neutron 
gas could exist as a stable configuration in neutron stars 
where they are bounded by gravitation, while isolated 
nuclei that exist for instance on earth are limited to the 
drip lines. 

The first prediction of the suppression of pairing in 
overflowing Z ^ 50 nuclear systems was performed in 
Ref. [^. It was proposed to attribute this suppression 
to the large coherence length of the weakly-superfluid 
neutrons gas: the neutron gas can penetrate the dense 
nuclear system and could impose its weak pairing field. 
It was however also noted in the same work that neutrons 
are dripping out of the double-magic nucleus ^^^Sn. In 
a recent work, the pairing gap of the last occupied state 
was also predicted to be quenched in overflowing nuclear 
systems, metallic grains and cold- atoms [5]. It was con- 
cluded that the suppression of superfluidity is a generic 
fact of a fermionic superfluid overflowing from a narrow 
container into a much wider one. It shall be noted that 
in this work the nuclear model was mainly limited to 
£=0 single particle states (s-states). The description of 
nuclear systems requires however the inclusion of ^ > 



single particle states for the bound and the mean-field 
resonance states. It is indeed well known that, close to 
the drip lines, resonance states play an important role 
for pairing properties and it is referred generally as con- 
tinuum coupling |8|, |9j. It has recently regained some 
attention in nuclei close to the drip-line (see for instance 
Refs. [iO: [ill and references therein). The general ques- 
tion of continuum coupling in overflowing superfluid sys- 
tems and various phenomenons such as persistence, sup- 
pression and reentrance of pairing remains to be stud- 
ied. Pairing correlations in the ground state of weakly- 
bound nuclei are commonly described by the Hartree- 
Fock-Bogohubov (HFB) theory [H, [ll. In most of the 
HFB calculations the continuum is discretizcd by solving 
the HFB equations with box boundary conditions jlj] . In 
this Letter a systematic analysis based on several over- 
flowing isotopic chains is performed and it is shown that 
pairing quenching and continuum coupling are strongly 
related. A pairing reentrance phenomenon with increas- 
ing temperature is predicted for Z = 50 overflowing sys- 
tems. 

In the present work, an HFB approach in coordinate 
representation is employed. This model has already been 
applied to describe nuclei and Wigner-Seitz cells in a 
fully self-consistent framework (see Ref. [ia| and refer- 
ences therein). The Skyrme SLy4 interaction 16] is used 
in the mean-field channel, and is completed by the ISS 
pairing force which is adjusted to the BCS pairing gap 
predicted by bare nucleon-nucleon potentials [la]- All 
the bound states are considered, the angular momentum 
goes up to Jmax — 27/2, and we take a large box ra- 
dius Rbox—'2'6-7 fm ensuring an adequate description of 
the continuum. An almost constant neutron density at 
the edge of the Wigner-Seitz cell is obtained using mixed 
Dirichlet-Von Neumann boundary conditions. The ex- 
pected nuclei in neutron stars located at the transition 
between the outer crust and the inner crust have a proton 
number Z around 30 to 50, depending on the models |l7| . 
We therefore selected 8 isotopes given in Table |T] and lo- 
cated around Z=28, 40 and 50. It is well adapted to 
perform a spherical HFB calculation since these nuclei 
are predicted spherical near the neutron drip-line [l8|. 

For these selected isotopes, the neutron-drip number 
Ndrip^ defined as the neutron number of the last nucleus 
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TABLE I: Isotope acronym, number of protons Z, number of 
neutrons of the last nucleus before the drip line Ndrip, for the 
selected set. Isotopes for which the drip-line nucleus is non- 
magic (magic) belongs to the group Ai {A2, respectively). 
The total occupation number of resonance states for the drip- 
line nuclei Nres is shown in the last column (see the text for 
more details). 
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FIG. 1: (color online) Neutron pairing gaps versus neutron 
density A'^ (top panel) and versus N~N(drip) (bottom panel) 
for the isotopes in the group ^1 (see the text for more details). 



before the two- neutrons separation energy S2n{N) — 
E{N) — E{N — 1) changes its sign [iJl, is given in ta- 
ble m The total neutron occupation number of resonance 
states Nres for the drip-line nuclei is also given. The cal- 
culation of Nres requires the identification of resonance 
states which, in our case, are defined as positive energy 
states with rms radii lower than 10 fm. It is interesting 
to remarck from Table U that the drip-line nuclei with a 
number of neutrons that coincides with the well-known 
magic numbers (82 and 126) have no states in the con- 
tinuum. The shell occupation in this case is or 1, as in 
magic nuclei. We therefore group Kr, Sr, Sn and Te iso- 
topes in the same group A2 ■ The persistence of magicity 
with the same magic numbers as in stable nuclei might 
be due to the sphericity of these nuclei. In the case of Ni, 
Zr, Mo and Ru, the drip-line occurs in more complicated 
shell structure partially involving resonance states and 
these isotopes are grouped in Ai. 

The drip-line isotopes identified in Table U can be con- 
sidered as the seed nuclei from which the overflow occurs. 
The dripping of neutrons shall therefore be influenced by 
the microscopic structure of these seed nuclei. Nuclei 
belonging to the group Ai have a resonance occupation 
number Nres which goes from 2 to 8 particles. This num- 
ber reveals the large continuum coupling. The number of 
particles in the resonance states is zero for nuclei belong- 
ing to the group ./I2, since in the magic nuclei located 
at the drip-line, pairing correlations are quenched and 
continuum coupling therefore largely suppressed. From 
a microscopic analysis, the structure of drip-line nuclei is 
qualitatively different between the group Ai and A2- In 
the following, we show how these differences impact the 
quenching of pairing correlations at overflowing. 

Figs. [Hand [2] display the average neutron pairing gaps 
versus the neutron number (top panel) and the differ- 
ence between the neutron number and the neutron drip 



N — Ndrip (bottom panel). The average neutron pairing 
gap is obtained from the neutron pairing fleld given by 
the HFB solution. Fig. [T] represents the neutron pair- 
ing gaps for the Ai isotopes while in Fig. [2] are shown 
the neutron pairing gaps for the A2 isotopes. In the 
case of nuclei from the Ai group, the neutron pairing 
gap persists beyond the drip line (see bottom panel of 
Fig. [1]). The continuum coupling preserves the pairing 
diffusivity around the Fermi energy: the occupancy of 
the scattering states beyond the drip-line does not sup- 
press superfuidity. The presence of an overflowing neu- 
tron gas has therefore a limited effect on the pairing gap 
for the Ai isotopes, showing the persistence of superflu- 
idity. In the case of the A2 isotopes, a suppression of the 
average pairing gap just beyond the drip Une is observed 
(Fig. [2]). The presence of a magic nucleus at the drip line 
have therefore a strong influence on the pairing gap at 
and beyond the drip line: the large spacing between the 
last occupied state and the first excited state (resonance 
state) suppresses superfiuidity. Figs. [1] and [2] provide a 
clear illustration that the shell structure of the drip-line 
nucleus determines the suppression or the persistence of 
superfluidity, even in overflowing nuclear systems. 

In Ref. [5|, the quenching of pairing in overflowing 
systems has been predicted upon overflow of trapped 
fermions. It should be underlined that in the present 
work, the mean-field in nuclear systems have a centrifu- 
gal term which gives rise to the mean-field resonance 
states and continuum coupling, whereas the conclusions 
obtained in Ref. [5] are applicable only to overflowing 
systems involving £ = single particle states. In more 
general systems, like in the crust of neutron stars, a more 
sophisticated theoretical model such as the present HFB 
one predicts a suppression of superfluidity only if the con- 
tinuum coupling is quenched, by the shell structure for 
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FIG. 2: (color online) Same as Fig. [T] for isotopes in the A2 
group. 
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FIG. 3: (Color online) Temperature averaged neutron pairing 
gap versus temperature for I60,i76,i80,200gj^^ ^ superfluid re- 
entrant effect is observed for ^^^Sn and ^*°Sn. 



instance. 

The persistence of superfluidity upon overflow of 
bound neutrons might not be the only consequence of 
continuum coupling: in the case where pairing is sup- 
pressed, the increase of temperature may generate the 
reentrance of superfluidity. The finite-temperature HFB 
model [231 , is employed to study the reentrance of pair- 
ing in the thermal state of overflowing even-even nu- 
clear systems. The temperature-averaged pairing gap for 
I60,i76,i80,200sn is showu in Fig.H In i^OSn and ^ooSn it 
behaves as expected from HFB theory: the pairing gap 
vanishes at the critical temperature Tc = 0.57A(r « 0) 
(see Ref. [2^ and references therein). In the case of 
^^^Sn and ^^''Sn, the reentrance of superfluidity is ob- 
served with increasing temperature. This reentrance is 
induced by the presence of resonances states in the spec- 
trum of these nuclear systems: Being slightly too high 
in energy, these states are not occupied at zero tempera- 
ture (see Table [H, while at finite temperature, they can 
be partially occupied from the Fermi-Dirac distribution. 
At low temperature, the pairing correlations can there- 
fore be switched on allowing the reapparition of the su- 
perfluid state. The reentrance critical temperature de- 
pends on the step in energy between the last occupied 
bound state and the first resonance one, which changes 
from one system to another, as observed in Fig. [3] for 
^^^Sn and ^®°Sn. In the case of ^™Sn, this energy step 
is too large to give rise to the reentrance of superflu- 
idity before the highest critical temperature is reached. 
The ^^'^Sn overflowing system has an interesting phase 
diagram including three critical temperatures: with in- 
creasing temperature, two of them correspond to the van- 
ishing of superfluidity (Tci ~ H keV and T^s ~ 1 MeV) 
and one to its reappearance {Tc2 ~ 300 keV). The lowest 
critical temperature Td is associated to the transition 
from the superfluid to the normal state in the overflow- 



ing neutron gas. The highest critical temperature Tcs is 
similar for i6oa76,i80gj^^ indicating that superfluidity has 
been restored in the seed nucleus of i^e.isogj^ between 
Tc2 and Tc^. This superfluidity is mainly built on reso- 
nances populated at finite temperature. More generally, 
pairing reentrance in hot systems is observed for isotopes 
belonging to the group A2 where resonances are too high 
in energy to participate to pairing at zero temperature 
but close enough to the last occupied state to be reached 
at finite temperature. 

Reentrance of superfluidity at finite temperature have 
been predicted in nuclear systems such as in odd- 
nuclei :19i], rotational motion of nuclei [20|, and the 
deuteron pairing channel in asymmetric infinite mat- 
ter [21[. It was aslo predicted in polarised '^''*He |22| - 
|24| and in spin asymmetric cold atom gas [23, [la]. In 
all these systems, pairing at zero temperature is gener- 
ated by an attraction among Fermions of different spin 
or isospin. Superfluidity is therefore maximum in spin 
or isospin symmetric systems for which there is a match- 
ing of the Fermi levels of the constituent Cooper-pairs. 
Breaking the spin or isospin symmetry disfavor pairing 
while temperature in asymmetric systems acts in favor of 
restoring the broken symmetry and can eventually induce 
a reentrance of pairing. At variance with this mechanism 
the pairing reentrance phenomenon discussed in this Let- 
ter is based on a novel mechanism in finite systems where 
resonance states play a major role. 

In summary, we have investigated the pairing proper- 
ties of nuclear systems upon overflowing superfluid neu- 
trons. Suppression, persistence and reentrance of super- 
fluidity can occur in these finite systems. From a sys- 
tematic HFB calculations on 8 isotopic chains, the pair- 
ing properties is shown to be strongly correlated to the 
continuum coupling, both at zero and finite temperature. 
At zero temperature, the coupling between the seed nu- 
cleus and the gas is weak, and a formal separation of their 



properties into a nucleus plus a gas provides a qualitative 
understanding of the suppression and the persistence of 
superfluidity. With increasing temperature in the nor- 
mal state, the Fermi-Dirac distribution can populate the 
resonance states giving rise to the reentrance of superflu- 
idity. The pairing correlations in the nuclear system are 
switched on again and consecutive critical temperatures 
are predicted. 

The understanding of the suppression, persistence and 
reentrance of superfluidity in nuclear systems, deeply re- 
lated to the continuum coupling, opens wide perspectives 
for discoveries in weakly bound nuclei, as well as it sheds 
new light on the transition between the inner and outer 
crusts in neutron stars. The role of resonances around 
the neutron drip not only changes the microscopic un- 
derstanding of the neutron drip-out mechanism, but it 
also modifies the thermal properties of the crust through 
the strength of the pairing interaction [23] . The temper- 
atures at work during cooling are typically of the order of 



10 to 500 keV [V\ and coincide with the critical tempera- 
tures of the reentrance phenomenon. As a consequence, 
the novel pairing reentrance phenomenon analyzed in this 
Letter is expected to modify the thermodynamical and 
cooling properties in the crust of neutron stars. The links 
with other superfluid Fermi-systems shall be investigated 
in the near future. For instance in cold Fermionic atoms 
overflowing from an inner trap to a larger one [3, |3| , it 
is known that stable path close to the centrifugal energy 
are the classical analog of the quantal resonances. It will 
therefore be interesting to investigate the role of these 
stable paths on the superfluid properties of overflowing 
cold Fermionic atoms. 
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